1 

Flocking of Multi-agent Dynamical Systems 
Based on Pseudo-leader Mechanism 

Jin Zhou, Wenwu Yu, Xiaoqun Wu, Michael Small, Jun-an Lu 

Abstract 

Flocking behavior of multiple agents can be widely observed in nature such as schooling fish 
and flocking birds. Recent literature has proposed the possibility that flocking is possible even only 
a small fraction of agents are informed of the desired position and velocity. However, it is still a 
challenging problem to determine which agents should be informed or have the ability to detect the 
desired information. This paper aims to address this problem. By combining the ideas of virtual force 
and pseudo-leader mechanism, where a pseudo-leader represents an agent who can detect the desired 
information, we propose a scheme for choosing pseudo-leaders in a multi-agent group. The presented 
scheme can be applied to a multi-agent group even with an unconnected or switching neighbor graph. 
Experiments are given to show that the methods presented in this paper are of high accuracy and perform 
well. 
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I. Introduction 

Flocking is a collective behavior of a large group of mobile agents. Typical flocking phenomena 
include flocks of birds, schools of fish, herds of animals and colonies of bacteria. In nature, 
animals achieve flocking for various reasons. For example, in order to keep safety in numbers 
and also to confuse predators, they will form flocks for protection. An agent is more likely to 
be attacked if it strays away from the flocking group. 

As a common demonstration of emergence and emergent behavior, flocking was first simulated 
on a computer by Craig Reynolds [1]. In 1987, he started with a boid model to build a simulated 
flock and introduced three rules to simulate flocking: 

• Collision Avoidance: steer to avoid collision with nearby flockmates (short range repulsion). 

• Velocity Matching: steer to match velocity with nearby flockmates. 

• Flock Centering: steer to stay close to nearby flockmates (long range attraction). 

The mechanism, known as "separation", "alignment" and "cohesion", results in all agents moving 
in a formation with the same heading and a fixed network structure. From then on, these three 
rules have been widely used to study flocking behavior. 

As a special case of Reynolds' model, in [2], Vicsek et al. presented a simulation model based 
on nearest neighborhood law, in which each agent's heading is updated by the average of the 
headings of its nearest neighbors and itself. It was shown that the headings of all the group 
agents converge to a common value. 

A lot of works have been published based on Reynolds' and Vicsek's models in recent 
years [3]-[ll]. To thoroughly and systematically investigate flocking behavior, artificial potential 
functions (APFs) are widely used. In [4], Tanner et al. presented an APF in a network with fixed 
topology which is a differential, nonnegative and radially unbounded function of the distance 
between two agents. Then, he modified the APF to a nonsmooth one in a network with switching 
topology which captures the fact that there is no agent interaction beyond a proper distance [5]. 
Later, Olfati-Saber modified Tanner's APF and defined another APF which is a bounded and 
smooth one for switching topology. 

In view of the pitfall of regular fragmentation [7], which is a phenomenon of flocking failure 
most likely occurring for generic set of initial states and large number of agents, Olfati-Saber 
also introduced a flocking mechanism based on a virtual leader [7]. Even though the initial 
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states is selected randomly, the mechanism can guarantee a flocking behavior. However, this 
method assumes that all the agents know information about the virtual leader. In consideration 
of practical use, Shi et al. and Su et al. showed that flocking appears even when only some 
agents are informed [8]-[10], these we call pseudo-leaders in this paper. 

Though Shi et al. and Su et al. showed that flocking could appear when not all the agents are 
informed in a group, precisely which agents should be informed has not yet been considered. 
In this paper, we focus on investigating which agents should be selected as pseudo-leaders for 
flocking in a weighted network. Another difference from previous work is that, in the absence 
of control, the acceleration of each agent is dynamic. This can be seen in our system model 
in Section II. That is, if there are no attractive/repulsive control, no information exchange with 
others and no information received from the virtual leader for an agent, its acceleration is dynamic 
instead of constant, as is the case in previous work. By using Lyapunov stability theory, a simple 
criterion for choosing pseudo-leaders is proposed. 

The reminder of this paper is organized as follows. In Section II, model depiction, preliminaries 
about graph theory, and mathematical analysis are briefly introduced. In Section III, the main 
results are proposed. We investigate how to select pseudo-leaders for flocking in a multi-agent 
dynamical group with fixed topology. Section IV gives some extensions and discussion of 
the main results so that one can gain useful insight into the problem of choosing pseudo- 
leaders. Two computational examples in small-sized and large-scale groups are simulated to 
illustrate effectiveness of the proposed approach in Section V. We summarize the main ideas and 
conclusions in Section VI. 

II. Preliminaries 

A. Graph theory 

To make this paper self-contained, some basics of graph theory are recalled [12]. 

A graph Q is a pair of sets (V, £), where V is a finite non-empty set of elements called 
vertices, and £ is a set of unordered pairs of distinct vertices called edges. The set V and £ 
are the vertex set and edge set of Q, and are often denoted by V(G) and 6(G), respectively. If 
i,j G V and (i,j) G 6, then % and j are adjacent vertices, or neighbors. The set of neighbors 
of a vertex z is Mi. A walk in a graph is a sequence of vertices and edges j , e\, j±, • • • , e^, jk, 
in which each edge e r = (j r -i,jr)- A path is a walk in which no vertex is repeated. If there is 
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a path between any two vertices of a graph Q, then Q is said to be connected. If Q and Q 1 are 
graphs with V{Qi) C V(<?) and C £((?), then & is a subgraph of 6?. 

The position and velocity neighbor graph Q = (V, £, W) is a weighted graph consisting of a 
set of indexed vertices and a set of ordered edges, where W = (wy) G R NxN is the weight 
matrix which represents the weighted coupling coefficients of interaction between the agents. If 
there is a link from vertex % to vertex j (j ^ i), then Wij = Wji > and Wij is the weight; 
otherwise, = 0. Throughout the paper, assume that W is a symmetric matrix satisfying 
diffusive condition 

N 

i=i 

B. Model depiction 

Consider a multi-agent system consisting of N agents. Here, the moving model of each agent 
in the group is given by 

Yi = f(Vj) + Ui, 

where 1 < % < N, p i e M™, G IR n and f(vj) G IR n denote the position, velocity and the 
acceleration dynamics (without control input) of the i-th agent respectively, and G W n is the 
control input of agent i. 

The motion model for the virtual leader is 

*■ = (2) 
= f(vj)> 

where p ; G M. n , V; G W l and f(vj) G R™ represent the position, velocity and the acceleration 
dynamics of the virtual leader, respectively. 

The flocking task is to design an appropriate control input Uj such that Reynolds' rules are 
followed, and then all agents moving in a formation with a common heading and collision 
avoidance. 
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C. Mathematical Preliminaries 

In this subsection, some useful mathematical definitions, lemmas and assumptions are outlined. 

Definition 1. A matrix A = (a^) G R N x N is called reducible if the indices 1, 2, • • • ,N can be 
divided into two disjoint nonempty sets ■ ■ and ji,j 2 , • • • ,jp 2 (with (3\ + f3 2 = N) 

such that 

for 1 < ki < /3i and 1 < k 2 < /?2- A matrix A = (a^) G IR 7VxAr is irreducible if and only 
if it is not reducible. 

A matrix is reducible if and only if it can be placed into block upper-triangular form by 
simultaneous row/column permutations. In addition, a matrix is reducible if and only if its 
associated graph is not connected [14]. 

Definition 2 [13] . A matrix A = (a^) G M. NxN is called diagonally dominant if |ajj| > 

N ' N 

\a>ij\ for 1 < i < N. A is called strictly diagonally dominant if |ajj| > Yl \ a ij\ 
for 1 < % < N. 

The Gershgorin circle theorem [13] results in many interesting conclusions. A strictly di- 
agonally dominant matrix is nonsingular. A symmetric diagonally dominant real matrix with 
nonnegative diagonal entries is positive semi-definite. If a symmetric matrix is strictly diagonally 
dominant and all its diagonal elements are positive, then its eigenvalues are positive; if all its 
diagonal elements are negative, then its eigenvalues are negative. Thus it is obvious that the 
real part of eigenvalues of the weight matrix W which are all negative except an eigenvalue 
with multiplicity one. Throughout the paper, we denote a positive definite (positive semi-definite, 
negative definite, negative semi-definite) matrix A as A > (> 0, < 0, < 0). 

Definition 3 [15] . A set is said to be a positively invariant set with respect to an equation 
if a solution x(t) of this equation satisfies 

x(0) G D, =>- x(t) G D, Vt > 0. 
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Definition 4 [15] . A continuous function g : [0, jo) — > [0, oo) is said to belong to class /C if it 
is strictly increasing and g(0) = 0. It is said to belong to class K.^ if p = oo and g(p) — > oo as 
p — > oo. 

Definition 5. 2-norm, also called Euclidean norm, of a vector £ is defined as 



where er is a positive constant. 

Due to the fact that ||f || ff = as ||£|| 2 = 0, \\£\\ a ^ oo as ||£|| 2 -> oo, and that §|||||| = 

, = = > when ||£|| 2 ^ 0, ||£|| CT is a class fC^ function of ||£|| 2 . On the other hand, note 

y l + (J \\i\\2 

that even though ||£|| CT is not a norm in the sense of algebra, it is differentiable everywhere while 
|| £|| 2 is not. Thus \\i\\ a instead of ||£|| 2 is used to construct APF in this paper. 

Definition 7. Artificial Potential Function (APF) Vij is a differentiable, nonnegative, radially 
unbounded [15] function of ||p i — Pj\\ a , which is the cr-norm of the position error between the 
i-th and the j-th agents. V^(||pj — Pj\\ a ) has the following properties: 
(!) ^-(llPi-PjIU) ^ oo as HPi — Pj|U -> 0, 

(2) V^j(||pj — poller) attains its unique minimum when the 2-th agent and the j-th agent are 
located at a desired distance. 

An example of APF is 



uh = Vet 



where T denotes the transpose of a vector or a matrix. 



Definition 6 [7] . a-norm of a vector £ is defined as 




VijdlPi -Pj\\a) = ci In llPi-Pjl^ + 




(3) 
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where c x > 0, c 2 > are two constants. The potential function approaches infinity as ||p i — p-^ 
tends to 0, and attains its unique minimum when ||Pj — Pj|| CT = yj^- 

In order to derive the main results, the following lemmas and assumption are needed. 

Lemma 1. (Schur complement [16], [17]) The following linear matrix inequality (LMI) 

I A{x) B{x) \ < Q 

where A(x) T = A(x), C(x) T = C(x), is equivalent to either of the following conditions: 

(a) A(x) < and C(x) - B(x) T A{x)~ l B{x) < 

(b) C(x) < and A(x) - B(x) C(x)" 1 B(x) T < 0. 

Lemma 2 [18] . If a symmetric matrix A = (a^) e M JVxAr {a^ > 0, % ^ j) is irre- 
ducible and satisfies diffusive condition, then A — Dj < holds for a diagonal matrix Dj = 
diag{0, ■ ■ ■ , 0, 5, 0, • • • , 0} G R N x N , where the 2-th (1 < % < N) element 5 is any positive 
constant and the others are 0. 

Assumption 1 (Al). Suppose that there exists a positive constant a satisfying (£ 2 — £i) T (f(£2) — 
f(6)) < a ||^2 - for any two vectors f 1} f 2 e R n . 



III. Main results 

In this section, a group of mobile agents with fixed topology is considered. Assume that the 
position and velocity neighbor graph Q = (V, 5, W) is connected. That is, the weight matrix 
W is irreducible. For the case that Q = (V,£, W) is unconnected, discussion will be found in 
Section IV. C. 

Generally, there should be attractive and repulsive mechanisms in the control input, where the 
attraction indicates that each agent wants to be close to nearby agents and repulsion provides 
the fact that each agent does not want to be too close to nearby flockmates. These two mecha- 
nisms can be jointly embodied in APF (actually there are many ways to achieve attraction and 
repulsion). 

Besides "separation" and "cohesion", "alignment" is an important rule in flocking. We use 
the approach that group members receive moving information from a virtual leader to realize 
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alignment. Based on this approach, recent literatures have shown the phenomenon that flocking 
will appear even if not all the agents are informed in the group. Here, we investigate which agents 
should be selected as informed ones for flocking in a weighted network. An agent, which utilizes 
the motion information of the virtual leader as a reference in its controller, is a pseudo-leader 
of the group. Suppose that the ii-th, i 2 -th, ■ ■ ■ , i m -th agents are chosen as the pseudo-leaders. 
The control law for the i-th agent is given by 



Ui = u] + uf + uf, 

N 

= - E WijiVi-Vj) = E WijVj, 

jeMi 3=1 

N 

«f = - E «'ijVp.^(||p i -p j y = - E Wij Vp^idlP; -Pj\\a), 

, f -^v(vv-v z ), h ir = k ir (\\y ir -vj|||) % e {i r \l < r < m, 1 < i r < N}, 
I otherwise, 

(4) 

where V represents the gradient of a function, h ir is the adaptive feedback gain, k ir > (1 < 
r < m, 1 < i r < N) is a constant. For agent uj is the velocity coupling of interactions 
between agents, and stands for the coupling gradient of APF with respect to its position. The 
term uf, which will disappear if agent % is not a pseudo-leader, controls the received information 
from the virtual leader. It provides an adaptive feedback adjusting mechanism and will be more 
practical in engineering than those with linear feedback ones. 

Define the position and velocity error between agent % and the virtual leader as ef = Pj — p^, 
ej = Vj — V;, then we have 

e p = e v 



Since Al holds, we obtain 



ej = f( Vi ) - f(v,) + u, 



eI T (f( Vi )-f(vO) < aM\\l 



for 1 < i < N. In addition, it is easily to get 
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V p ^(||p;- Pj y = V^Vydlef-eJIU) 

and 

VeP^dlef-e^lU) = -V^-dlef-eflU), 

and then it follows 

I^Hp.-p.y = 2ef T V e P^(||ef- ej p || CT ) 
= 2eI T V Pj ^(||p i -p J .|U). 

Theorem 1. Suppose that Al holds. If Wn-th + aljv-m < 0, where WAr_ m is the minor matrix 
of the weight matrix W by removing all the z r -th (1 < r < m, 1 < i r < N) row-column 
pairs, flocking behavior appears in system (1) by the control strategy (2) and (4). That is, the 
velocities of all agents approach the desired velocity asymptotically, collisions between agents 
are avoided and the distances between all agents are invariant. Furthermore, the global potentials 
of all agents in the group are minimized with the final configuration. 

Proof. Consider a positive semi-definite function as 

N N N 

l = \ E E ^ ^(IIp.-p.-IU) + \ E ef el 

i i V ( hir ~ h > 

r=l 

where h > is a constant to be determined. 

We now prove that f2 c = {(P; — p 3 -, ej) | L < c, c > 0}, the sub-level set of L, is compact. 
Firstly, from L < cone gets ||ej||| = eJ T eJ < 2 c and w i3 ;Vy(||pj — Pj|| CT ) < 2 c fori < i,j < 
AT. Due to the fact that the continuous function Vij is radially unbounded, V^(||pj — p ■ \\ a ) — > oo 
as ||p i — Pj\\ a — > 0, and that ||p i — poller is a class function with respect to ||p- — p - 1| 2 , there 
exists a positive constant d such that Hp, — p-|| 2 < c'. Thus fi c is a bounded set. Secondly, 
because il' c = {L\ L < c, c>0}is closed, f2 c is a closed set for the continuity of function 
L. Then according to Heine-Borel Theorem [19], f2 c is compact. 
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Regarding L as a Lyapunov candidate, its derivative along the trajectories of (1), (2) and (4) 



is 



N I N N 

E eI T ffa) - !(▼,) + E e ) " E «; y V Pi ^(|| Pi - Pj \\ a ) 

N N mm 

+ E E ^ef V Pi ^(|| Pi - Pj .|| CT ) - E ^eT T eT + E - &K T eT 



< E aeT ej + E E ^ej 1 ej - E 1 e? r 

i=l i=l j=l r=l 



= E alkilll + E ^ll e IIII + E E Wij^Je] - e MKIII 

i=l j=l i=l j=lj^i r=l 

AT AT AT AT m 

< E « Il e lll2 + E ll e IIII + E E Mh \W)h - h \K\\i 

i=l i=l i=l j=ljy£i r=l 

= e vT Qe v , 

where e v = ( 1 1 ©a 1 1 2 , 1 1 1 1 2 , - - - , 1 1 e jv 1 1 2 ) , Q = olIn + W — H, and H is a diagonal matrix 
whose z r -th (1 < r < m) elements are h and the others are 0. 
After applying row-column permutation, Q can be changed into 

W* + al m - hl m W** 

W** T W N - m + aI N . m 

where W*,W** are the corresponding matrices with compatible dimensions. Since Wjv_ m + 
al N _ m < 0, W j\j — m + aI N ^ m is invertible. According to Lemma 1, choosing h be a positive 
constant satisfying W* + al m - hl m - W** (W N - m + a\ N - m )~ l W** T < 0, we have Q < 0. 
Furthermore, L < 0, L(t) is a non-increasing function of t. Thus any solution of (1), (2) and 
(4) starting in f2 c will stay in it. Namely, fi c is a positively invariant set. 

Because f2 c is compact and positively invariant, every solution of the system converges to the 
largest invariant set f2* of the set {(p i — Pj,eJ) | L = 0} on the basis of LaSalle's invariance 
principle [15]. In ef = ej = 0, which means that all agent velocities are equal and their 
position differences remain unchanged in steady state. 

Furthermore, ej = leads to ej = in Q*. Combining with equation (1), (2) and (4), we 

N 

have E w ij^Vi^ij(\\Pi ~ PjlU) = for 1 < i < N. This indicates that the group final 
configuration is a local minima of global potential function of agent i. 
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Collision avoidance can be proved by contradiction. Assume that there exists a time ti > 
so that the position difference between two distinct agents i* and i** satisfies ||pj* — Pj*» H2 — > 
as t — > ii. Then ||p it — Pj**|| CT — > since || . || CT is a class /C^ function with respect to || . || 2 . 
According to the definition of APF, V^(||pj* — Pj»»|| CT ) — > 00. This is in contradiction with the 
fact that tt c is a positively invariant set. Therefore, no two agents collide at any time t > 0. 

Thus the proof is completed. □ 

From this theorem, if X ma x(^N-m) < where X max (.) represents the maximum eigenvalue 
of a symmetric matrix, V = i 2 , • ■ ■ , i m } can be chosen as the pseudo-leader set to guarantee 
flocking in system (1), (2) and (4). 

IV. Discussion and Extensions 

In this section, some remarks and extensions are discussed to give some insights into the main 
results. 



A. Other flocking models 

Previous papers studied a classical model for flocking behavior in such as [4]-[9]. The motion 
models of the i-th agent and the virtual leader respectively are 

*' = V " (5) 

and 

* = V " (6) 
v/ = 0, 

where p i5 Vj (i = 1, 2, • • • , N, I), u«(i = 1, 2, • • • , N) and f have the same meaning with the 
equations (1) and (2). It is assumed that if there are no attractive/repulsive control, no information 
exchange with others and no information receiving from the virtual leader for an agent, then its 
acceleration is f(vj) = f(v/) = 0. For the virtual leader, f(v/) = implies that it moves along 
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a straight line with a desired velocity V;(0). Letting a = 0, we have the following theorem for 
choosing pseudo-leaders in the classical model. 

Theorem 2. If Wat-™ < holds, flocking behavior will appear in system (5) by the control 
strategy (4) and (6). That is, the velocities of all agents approach the desired velocity asymptoti- 
cally, collisions between agents are avoided and the distances between all agents are invariant. In 
addition, the global potentials of all these agents are minimized with the group final configuration. 

For other flocking models such as taking velocity damping into consideration [10] (which is 
frequently unavoidable when objects move with high speeds or in a viscous environment), the 
schemes for determining pseudo-leaders can be attained by similar analysis. 

B. A single pseudo-leader is enough for flocking 

Consider the classical model at first. Denote the minor matrix of W by deleting any row- 
column pair as Wjv-i- It can be rewritten as Wjv-i = Wjv_i + ^2~Di, where 1 < i < x-> 

i 

1 < X < N — 1, Wjv-i is the corresponding symmetric and diffusive matrix, and Dj is the 
diagonal matrix where the i-th element is negative and the others are 0. According to Lemma 
2, Wjv_i is negative definite. Therefore, flocking will occur with just one single pseudo-leader 
(any agent is available as an option) based on Theorem 2. 

Below we will discuss the model presented in Section II. B. Suppose that the weight matrix 
W = wB, where w > is the common weight coupling, B = (6^) e M JVxAr is the adjacent 
matrix with bij = bji = 1 if there is a link from agent i to agent j (j ^ i) and 6^ = 

N 

otherwise. Moreover, B satisfies diffusive condition hj — 0. For an agent i r (1 < r < 

3=1 

m, 1 < % r < N), if the minor matrix B N _i that is obtained by deleting the i r -th row-column 
pair of B satisfies A max (Bjv_i) < — ^ , it can be picked out as a pseudo-leader for flocking 
according to Theorem 1. Since B N ~i is negative definite by Gerschgorin theorem, it is concluded 
that flocking will be achieved with just one single pseudo-leader, which can be selected randomly 
from the vertex set V, provided that the common coupling weight w is large enough. 
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C. The position and velocity neighbor graph is unconnected 

In the main results, we assume that the position and velocity neighbor graph Q = (V,£, W) 
is connected. In reality, however, it is not always the case. Suppose that the graph consists 

of several connected subgraphs G j = (V j ,£ j ,W) (j G {1,2, •••}), where [j V j = V, 

j 

1J £ J = £, W- 7 G R Nj x Nj is the weight matrix of subgraph Qi which is symmetric, irreducible 

3 

and diffusive, and Nj = N. In subgraph Q j , pseudo-leader set V j can be picked out according 

3 

to Theorem 1 or Theorem 2. Put the pseudo-leaders in all the subgraphs together, the scheme 
for choosing pseudo-leader set V = (J V j of the whole unconnected group is obtained. 

3 

As discussed in the previous subsection, one single pseudo-leader in each connected subgroup 
is enough for flocking. If, however, there exists a connected subgroup in which no agent is 
informed, it is impossible for this subgroup to detect the desired moving information. Thus 
flocking failure, such as regular fragmentation, may occur under this circumstance. 

D. The position and velocity neighbor graph is switching 

The aforementioned results are on the basis of a multi-agent group whose topology is fixed. 
We can also consider a switching position and velocity neighbor graph Q(i) = (V,£(t),W(t)), 
where 

W(t) = W(r fc ) t G [r k , r k+1 ), k = 0, 1, • • • 

with To = is a switching matrix. It is easy to see that the graph switches at some instant time. 
Since W(rfc) is invariant in [r k , r fc+i ), a temporal pseudo-leader set V(r fc ) in this time interval can 
be determined based on Theorem 1 or Theorem 2. Thus for a multi-agent group whose topology 
is switching, a switching scheme can be employed to choose pseudo-leaders for flocking. 

For those switching models which utilize local information (nearest neighbor law) to establish 
the group topology [7]-[9], this technique will be more efficient than using a fixed pseudo-leader 
set. 

E. Center of the group members 

Define the position and velocity of the center of all group members as 
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N 

- A 1 

P = N Is Pi' 

i=l 
AT 

- A 1 

V = AT E Vi. 

i=l 



Then we have the error system 



eP = p 



— V A — 



AT 

Pi = jv E e i ' 

i=l 

1=1 



According to the proof of Theorem 1, every solution of the system converges to the largest 
invariant set Q* of the set {(p i — Pj,e]) \ L — 0}, in which ef = e] = for 1 < i < N. 
This results in e p = e v = in Q*. That is, in steady state, the center of all agent velocities 
is equal to the desired velocity, and the position difference between the center and the virtual 
leader remains unchanged. 



V. Numerical Simulations 

A. Flocking in a small-sized multi-agent group 

In the following, Theorem 1 is illustrated by using Lii system [20] as the dynamical ac- 
celeration in system (1) and (2). As a typical benchmark chaotic system, Lii system is given 
by 



/ 



x = 



\ 



J 



%2 

V X3 J 



( 



+ 



\ 



o 

-x ± x 3 
x x x 2 



\ 



J 



—a a 
c 
0-6 
^ Rx + T(x) , 

which has a chaotic attractor when a = 36, 6 = 3, c = 20. For any two state vectors y and z 
of Lii system, there exist constants M s such that ||y s ||, ||z s || < M s for 1 < s < 3 since the 
Lii attractor is bounded within a certain region. From simple numerical calculation, Mi = 25, 
M 2 = 30, M 3 = 45 is obtained. Therefore, one has 
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(y-z) T (f(y)-f(z)) 



(y 



R(y 



(y-z) T (T(y)-T(z)) 



( ( -a § ^ 



= (y - z) T R (y - z) + (y 1 - z x ) (z 2 y 3 - y 2 z 3 ) 

= (y-z) T 



< A 



—a 



a 
2 





c 








+ 
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2 
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(y-z) 



ll(y-z) 



w 60.3402 ||(y-z)|| 2 
Thus Lii system satisfies the assumption Al with a = 60.3402. 
For simplicity, consider the group with ten nodes whose adjacent matrix is 
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The common weight coupling is w = 70. By deleting the first three row-column pairs of B, the 
maximum eigenvalue of the minor matrix is A maa; (B 7 ) = —1. It is easy to see that 



A, 



(B 7 ) < - 



a 



w 



Thus the first three agents can be picked out as pseudo-leaders in the group according to Theorem 
1. 
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Fig. 1. Initial moving state of the group when N = 10 and m = 3. 



Simulation results are shown in Fig. 1 - Fig. 4 . In the simulation, the initial positions of 
the ten agents in the group are distributed randomly from the cube [0, 5] 3 . The initial velocity 
coordinates are randomly chosen from the cube [0,2] 3 . The initial position and velocity of the 
virtual leader, which is marked with a red star in Fig. 1 and Fig. 2, are set as p z (0) = (6, 6, 6) T 
and \i(0) = (1,1, 1) T respectively. Other parameters are chosen as h ir (0) = 1 and k ir = 
0.01 (1 < r < m, 1 < % T < N). Letting a = 0.1 in a-norm and c\ = c 2 = |, we use the 
following APF 

Vij(\\Pi-Pjh) = ci In llPi-PjH' + I, ^ ||2 . 

The group's moving states at t = and t = 30 are illustrated in Fig. 1 and Fig. 2. Here, the 
solid circles represent the pseudo-leaders that receive moving information of the virtual leader, 
and the hollow ones denote the followers in the group. The arrows display velocity vectors of 
all the agents. The dash lines depict the connections between agents. From these two figures, it 
is seen that in spite of the initial disordered state, all the agents flock at t = 30. 

The state of position and velocity differences between the agents and the virtual leader are 
illustrated in Fig. 3 and Fig. 4 respectively. It is shown that the velocities of all the group 
members converge to the desired velocity. Moreover, the position errors between agents remain 
fixed after a period of time. 
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Fig. 2. Group's moving state when N = 10 and m = 3 at t = 30. 
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Fig. 3. The position differences e? s (1 < i < 10, 1 < s < 3) between the agents and the virtual leader when N = 10 and 
m — 3. 



B. Flocking in a 200-sized multi-agent group 

Large-scale groups with complex topology are common in nature. It has been demonstrated 
that the topological information of most large-sized systems display scale-free features, among 
which Barabasi and Albert (BA) model [21] of preferential attachment has become the standard 
mechanism to explain the emergence of scale-free networks. Nodes are added to the network 
with a preferential bias toward attachment to nodes with already high degree. This naturally 
gives rise to hubs with a degree distribution following a power-law. 

In this subsection, a BA scale-free network consisting of 200 agents with q = 5 and q = 5 are 
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Fig. 4. The velocity differences e\ a (1 < i < 10, 1 < s < 3) between the agents and the virtual leader when iV = 10 and 

m = 3. 



considered, where q Q is the size of the initial network, and q is the number of edges added in each 
step. Similar to the previous simulation, we take the Lii system as the dynamical acceleration 
in system (1) and (2). The initial positions and the velocities of the 200 agents are selected 
randomly from the cube [0, 50] 3 and [0, 2] 3 , and that of the virtual leader are selected as p^O) = 
(60, 60, 60) T and V;(0) = (1, 1, 1) T . Assume that the other parameters in the simulation are the 
same as those in the previous subsection. 

Denote as the minor matrix of the adjacent matrix B BA by removing 6 row-column 
pairs which corresponds to the 6 agents with the largest degree in the whole group. Since the 
maximum eigenvalue of B^J satisfies 

XmaAKi) = -0.9284 < -0.8620 = ~, 

the 6 agents with the largest degrees can be picked out as pseudo-leaders of the group. After 
using the control mechanism presented in Theorem 1, flocking appears by letting just 3% agents 
be informed. The moving states of all the group members at t — and t = 30 are exhibited in 
Fig. 5 and Fig. 6. 

Clearly, the approaches presented in this paper are of high accuracy with good performance 
not only for small-sized multi-agent groups but also for larger-scale multi-agent systems. 
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Fig. 5. Initial moving state of the group when N = 200 and m = 6. 




Fig. 6. Group's moving state when N = 200 and m = 6 at t = 30. 



VI. Conclusions 

In this paper, we have presented a criterion for choosing pseudo-leaders in a multi-agent 
dynamical group. Particularly, the weight configuration of the position and velocity neighbor 
graph is not necessarily irreducible or time invariant. By combining the ideas of virtual force and 
pseudo-leader mechanism, mathematical analysis has been deduced to illustrate how to determine 
the pseudo-leader set in a group. The proposed schemes have been proved rigorously by using 
Schur complement and Lyapunov stability theory. Finally, two computational examples including 
small-sized and larger-scale multi-agent groups have been shown to illustrate the effectiveness 
of the proposed approach. 
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